Chapter 2 
Classical Models 



In this Chapter we will discuss models that can be described by classical statistical 

mechanics. We will concentrate on the classical spin models which are used not 
only to study magnetism but are valid also for other systems like binary alloys or 
lattice gases. 

In principle, all particles obey quantum statistical mechanics. Only when the 
temperature is so high and the density is so low that the average separation be- 
tween the particles is much larger than the thermal wavelength A, 

(-) »A,whereA=(^— , (2.1) 

the quantum statistical mechanics reduces to classical statistical mechanics and we 
can use a classical, Maxwell-Boltzmann distribution function (instead of Fermi- 
Dirac or Bose-Einstein). 

For magnetic materials, the situation is somewhat different. The Hamiltonian 
of a magnetic system is a function of spin operators which can usually not be 
directly approximated by classical vectors. Quantum models in d dimensions can 
be mapped onto classical models in d' — {d + 1) dimensions. The classical limit 
of the Heisenberg model, however, can be constructed for large eigenvalues of 
the spin operator by replacing the spin operators by three-dimensional classical 
vectors. 

There is another quantum-mechanical effect we must discuss. In quantum me- 
chanics identical particles are indistinguishable from each other whereas in classi- 
cal mechanics they are distinguishable. Therefore, the partition function of a sys- 
tem of non-interacting, non-localized particles (ideal gas!) is not just the product 
of single-particle partition functions, as one would expect from classical statistical 
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mechanics, but we must take into account that the particles are indistinguishable. 
We must divide the total partition function by the number of permutations between 
N identical particles, N\, 

ZM^^{Z,f, (2.2) 

where Zi is the one-particle partition function. The factor A^! is a purely quantum 
effect and could not be obtained from classical statistical mechanics. 

In (insulating) magnetic systems we deal with localized particles (spins) and 
permutations among particles are not possible. Therefore we must not divide the 
partition function hy {N\). 

A Hamiltonian of a simple magnetic system is: 

H — —-2^ JijSiSj — H 2_^Si. (2.3) 

The first term is the interaction energy between the spins, Jij being the (ex- 
change) interaction energy between the spins at the sites i and j. The sum is 

over all bonds between the spins. S is the spin operator; in general it is a vec- 
tor. The order parameter is equal to the thermal and quantum-mechanical av- 
erage of S, m = {{S)). The thermal average of the first term in Eq. (2.3) is 
E{S, m) — {E{S, {'5'})). (Here, S is the entropy and {S} is a spin configuration. 
The thermal average has to be taken over all spin configurations.) E{S, m) is the 
internal energy of a magnetic system. 

The second term in (2.3) is the interaction with external magnetic field H. The 
thermal average of this term corresponds to PV in fluids. We see that {H{S, H)) 
is a Legendre transformation of -^(^S, m) and is thus equivalent to the enthalpy of 
fluids. 



2.1 Real Gases 

In the ideal gas we neglected the collisions between the particles. In reality, the 
atoms (molecules) have a finite diameter and when they come close together, they 
interact. The Hamiltonian of interacting particles in a box is: 



2m 



+ u{n) 



(2.4) 
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u{ri — r-j) = u{r) is the interaction energy and U the potential of the box walls. 
The particles are polarizable and this leads to an attractive van der Waals interac- 
tion between the induced electric dipole moments. Such an interaction energy is 
proportional to r~^, where r is the interatomic distance. At short distances, how- 
ever, the electrons repel because of the Pauli exclusion principle. Although the 
repulsive potential energy usually decays exponentially with r, we will write it in 
the form: 



u{r) 



oo 



-e{2ro/rY 



r < 2ro 
r > 2ro 



(2.5) 



(2ro is the minimal distance to which the particles of "hard-core" radius ro can 
approach.) 

The partition function of N interacting identical particles is 




Figure 2.1: Interatomic potential energy ^(r) in real gases. 
Z{T,V,N)^j^l.-- 1 d^p,---d^pN /•••/ o?Vi---ciVjvX 

After integrating over all momenta, the partition function becomes: 

Z{T,V,N) = ^(^'^)''''' Z,{T,V,Nl (2J) 

where Zj carries all the interactions between the particles. We calculate the par- 
tition function in a mean-field approximation in which all the correlations are ne- 
glected. The potential acting on one particle is calculated by assuming a uniform 
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spatial distribution of other particles. 



fVi „ 

/ d'^rie" 

^0 



N 



(2.8) 



Vi is the volume, available to one particle, Vi = V - (iV/2)(47r/3)(2ro)^ = 
V — Nb. Here we subtracted from V the volume occupied by other particles. N/2 
is the number of pairs of particles. The mean-field potential $ is: 



roo 

$ = - / e {2ro/rf{N/V)ATrr^dr = 



2Na 



(2.9) 



The resulting partition function is 

Z{T, V, N) = ^('^)'"" (V - iV6) V"""/" 



(a = (27r/3)e {2rof, b = (27r/3)(2ro)^), the Helmholtz free energy is: 
F(T, V, N) = -fcsTln Z = -NksT In 



V-Nb f2nmkBTy^'^ 



and the equation of state is 



p{T,V,N)^ 



NknT N^a 



'dV'T,N V-Nb V^' 
The last equation can be written in a more familiar form: 



(2.10) 



V ' 
(2.11) 



(2.12) 



N'^a 

{P+^){V -Nb) = NkBT. 



(2.13) 



This is the familiar van der Waals equation. We derived the van der Waals equa- 
tion of state in a mean-field approximation, we neglected all the correlations be- 
tween the molecules. Nevertheless, this equation qualitatively correctly describes 
the liquid-gas transition and it also predicts the existence of the critical point. 
However, as we neglected the correlations, which are essential in the vicinity of 
the critical point, it does not predict the correct critical behaviour. 
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The critical point is located at the inflexion point of the isotherm, {dp/dV)j. — 
{d'^p/dV'^)j: — 0, and is located at: 

Vc = ^Nh 
a 

" 2762 
8a 

= mk; 

Below the critical point, the liquid and the gas phases are separated by a first- 
order phase transition. The equilibrium transition between the two phases takes 
place when their chemical potentials and therefore the Gibbs free energies are 
equal. 

Some comments are in place here. 

• Critical exponents. The van der Waals equation was derived in a mean-field 
approximation, therefore all the exponents are mean-field like. 

• Liquid gas transition. The liquid phase is stable until the point A is 
reached. In equilibrium, at this point, liquid starts to evaporate, the system 
enters the coexistence region (mixture of gas and liquid) until the point B is 
reached, which corresponds to the pure gas phase. 

If equilibrium is not reached, the liquid phase persists beyond the point A 
towards the point D. In this region, the liquid is metastable and can (at 
least in principle) persist until the point D is reached. Beyond this point, 
the compressibility is negative, the liquid state is unstable (not accessible). 

• Gas liquid transition. The situation is reversed. In equilibrium, the 
gas phase is stable until the point B is reached. This is the onset of the 
coexistence region which exists until the point A is reached. Metastable gas 
phase exists between the points B and E and, again, the region E to I? is 
unstable. 



• Convexity. As we have discussed in Section 2.3.3 already, the Helmholtz 
free energy F of a stable system has to be a concave function of V. For 
metastable states it is concave locally whereas for unstable states, it is con- 
vex. 
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• Spinodals. The transition from a metastable (homogeneous) state toward 
the equihbrium (mixed) state is connected with the nucleation of droplets of 
the new phase. Droplets first nucleate and then grow until the equilibrium 
(coexistence) is reached between the new and the old phase. Nucleation is 
a dynamic process. 

The ultimate stabiUty Umit of the metastable state is called the spinodal Une. 
Along this line the isothermal compressibility diverges, therefore the spin- 
odal line is in a sense a line of critical points. In the van der Waals (mean- 
field) systems the separation between the metastable and unstable states is 
sharp. Beyond MFA, for systems with short-range correlations, however, 
the transition between metastable and unstable states is often gradual. 

2.2 Ising Model 

The Ising model is the prototype model for all magnetic phase transitions and is 
probably the most studied model of statistical physics. In this model, the spin 

operator S is replaced by a number, which represents the z— component of the 
spin and is usually 5 = ±1 ("up" or "down"). The order parameter m = (S) is 
thus one-dimensional, it is a two-component scalar. 
The Ising Hamiltonian is 

n = -W JijSiSj -HY^Si (2.15) 

Jij is the interaction energy between the spins i and j. Usually (not always) the 
sum is only over nearest neighbouring (NN) pairs, then — J. H is the external 
magnetic field (in energy units). The case J > corresponds to ferromagnetism 
and J < to antiferromagnetism. As we shall show later, the Ising model has a 
phase transition at finite T if d > 1. In li = 1 the transition is at T = 0. Therefore 
we say that d — lis the lower critical dimension of the Ising model. 

In the low-temperature phase, the symmetry is broken, the spins are ordered, 
the order parameter is m 7^ 0, whereas in the high-temperature phase the sym- 
metry is not broken, m vanishes, the spins are disordered. Therefore this model 
is a prototype of all order-disorder transitions. Other examples of order-disorder 
transitions, I would like to mention, are lattice gases and binary alloys. All these 
systems are equivalent to each other, they all have the same critical behaviour 
(which - of course - depends on d), for given d they all belong to the same univer- 
sality class. 
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2.2.1 Lattice Gas 

Imagine that a fluid (gas or hquid) system is divided into a regular lattice of cells 
of volume roughly equal to the particle (atom, molecule) volume. We say that 
the cell is occupied if a particle centre falls into this cell. Since the cell volume is 
comparable to the particle volume, there can be not more than one particle per cell. 
Notice that microscopic details (like positions of the particles) are irrelevant in the 
critical region, therefore the lattice gas model is suitable for studying the critical 
behaviour of real gases although the space is discretized. (It is also assumed that 
the kinetic energy of the gas particles, which is neglected in the lattice gas model, 
does not contribute a singular part to the partition function.) 

o o o o o 
o o o o • o o 
o •-• o o o o 
o i-i o o • o 
o o o o o • o 

Figure 2.2: A configuration of a lattice gas md — 2. Solid circles represent atoms 
and open circles empty sites. Lines show bonds between nearest-neighbouring 
atoms. 

A lattice gas is thus a collection of particles whose kinetic energy is neglected 
and that are arranged into discrete cells. The cells are either occupied by one atom 
[closed circles in Fig. 2.2] or empty [open circles]. The Hamiltonian is: 

where u is the interaction potential, 

J- oo {vij = 0) 
u{rij) = I -u (rjj = NN distance) (2.17) 
(otherwise) 

and {i.j) is the sum over pairs, pi is the occupation number, it is 1 if the site is 
occupied by an atom, and zero otherwise. The occupation is, like the Ising spin. 
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a double- valued function. Therefore we expect that the two models are related. 
Indeed, if we set 



the Hamiltonian 2.16 becomes: 



Here, z is the number of NN neighbours (4 for square lattice, 6 for simple cubic 
lattice), and TYo is a constant, independent of S. Thus, we see that 1-iig is identical 
to the Ising Hamiltonian if we set ii/4 = J and {ji/2 -\- uz/A) — H. For w > 0, 
the atoms attract and they "condense" at low temperatures whereas for m < 
the atoms repel. At 1/2 coverage (one half of the lattice sites is occupied) and 
M < 0, the empty and occupied sites will alternate at low T, like spins in an 
antiferromagnet. 

The lattice gas model in d = 2 is used to investigate, e.g., hydrogen adsorbed 
on metal surfaces. The substrate provides discrete lattice sites on which hydrogen 
can be adsorbed. By varying the pressure, the coverage (amount of adsorbed 
hydrogen) is varied. This corresponds to varying the field in magnetic systems. 

There is an important difference between the lattice gas and the Ising models 
of magnetism. The number of atoms is constant (provided no atoms evaporate 
from or condense on the lattice), we say that the order parameter in the lattice gas 
models is conserved whereas it is not conserved in magnetic systems. This has 
important consequences for dynamics. 

2.2.2 Binary Alloys 

A binary alloy is composed of two metals, say A and B, arranged on a lattice. We 
introduce the occupations: 




(2.19) 




1 




if i is occupied by an A atom 
if i is occupied by a B atom 




1 




if i is occupied by a B atom 
if i is occupied by an A atom 



Pi + qi = 1. 



(2.20) 
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Depending on which atoms are nearest neighbours, there are three different NN 
interaction energies: uaa, ubb, and uab — uba- The Hamiltonian is: 

T-Cba = -Uaa J2 PiPj ~ '^BB QiQj - Uab J2 (Pi^j + ^iPj)^ (^-^l) 

ihj) (iJ) 

where each sum runs over all NN pairs. If we write pi — {1 + Si)/2 and 
Qi — {1 — Si)/2 then Si — 1 if i occupied by A and Si — —1 for a B site. 
When there are 50% A atoms and 50% B atoms, Si = 0, and the Hamiltonian 
becomes an Ising Hamiltonian in the absence of a magnetic field: 

'Hba — ~J'Y SiSj + 7io- (2.22) 

(id) 

Ho is a constant, and J = {uaa + ubb — '^uab)/^- For binary alloys, the order 
parameter is also conserved. 

2.2.3 Ising model md=l 

It is not difficult to solve the Ising model in the absence of H m d = 1 exactly. 
The Hamiltonian of a linear chain of N spins with NN interactions is 

T-iid — SiSi+i. (2.23) 

i 

We shall use periodic boundary conditions, that means that the spins will be 




Figure 2.3: Lattice sites of an Ising model in one dimension with periodic bound- 
ary conditions. 
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arranged on a ring, see Fig. 2.3. Then, the A^— th and the first spins are NN and 
the system is periodic. The partition function is: 

TV 

Z = Y. e-^««^» = = E IT 6^-^^*^'+^ . (2.24) 

{S} {S} {S} i=l 

(Don't forget: J2{s} is the sum over all spin configurations - there are 2^ spin 
configurations - whereas J2i is the sum over all sitesl) 
The exponent in (2.24) can be written as: 

e/3JSiSi+i ^ cosh (3JSiSi+i+smh (3JSiSi+i = cosh (3J+SiSi+i sinh (3 J. (2.25) 

The last equality holds because SiSi+i can only be ±1 and cosh is an even function 
whereas sinh is an odd function of the argument. The partition function is: 

N 

Z = {cosh 13 J f E n [1 + > (2-26) 

{5} i=l 

where K = tanh (3 J. We work out the product and sort terms in powers of K: 

Z = {coshpjfJ2(^ + SiS2K){l + S2SsK) •••(! + S^SiK) 

{s} 

= {cosh 13 J f ^ [1 + K{SiS2 + S2S3 + --- + SnS,)+ 

{S} 

+K^{S,S2S2Ss + ...)+... + K''{SiS2S2 ■ ■ ■ SmSnS^)] . (2.27) 

The terms, linear in K contain products of two different (neighbouring) spins, Uke 
SiSi+i. The sum over all spin configurations of this product vanish, J2{s} SiSi+i = 
0, because there are two configurations with parallel spins (SiSi+i = +1) and two 
with antiparallel spins (SiSi+i — —1). Thus, the term linear in K vanishes after 
summation over all spin configurations. For the same reason also the sum over all 
spin configurations, which appear at the term proportional to K"^, vanish. In order 
for a term to be different from zero, all the spins in the product must appear twice 
(then, J2{Si} = 2). This condition is fulfilled only in the last term, which - af- 
ter summation over all spin configurations - gives 2^K^. Therefore the partition 
function of the Ising model of a linear chain of N spins is: 

Z ={2 cosh pjf \l + K^] . (2.28) 
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In the thermodynamic Umit, — > oo, — > 0, and the partition function be- 
comes: 



Z{H = 0) = (2cosh/5J) 



N 



(2.29) 



The free energy. 



F/N 



Temperature 



Figure 2.4: Free energy of the Ising model md — 1. 

F = -NkeT ln(2 cosh (3 J) , (2.30) 

is regular at all finite T, therefore the Ising model ind—l has no phase transition 
at any finite T. 

To investigate the possibility of a phase transition at T = 0, we will calculate 
the correlation function T and the correlation length ^. We first replace J by Jj 
(at the end we put Jj = J), so that 

Z = 2^[]cosh/3Jj-. (2.31) 

j 

The nearest-neighbour spin-spin correlation function is: 

^ {5} 

1 dZ 

The spin correlation function for arbitrary distance k, 

r{i,i + k) = (SiSi+k) 
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— {SiSi^iSi-^iSi+2Si-\-2Si-\-k-iSi+k-iSi+k) 
1 d d d d 

ZdiPJi) a(/?J,+i) d{(3Ji+2) ' ' ' d{(3Ji+k-i] 



Z, (2.33) 



V{k) = r(i, i + k) = [tanh (/5 J)]' 



(2.34) 



is just the product of all the nearest-neighbour correlation functions between the 
sites i and i + k, and is (if allJ's are equal) independent of the site i. The decay of 
the spin correlation function with distance is shown in Fig. 2.5. The correlations 
decay exponentially with distance and the correlation length is: 



1 



a In tanh /9 J' 



(2.35) 



where a is the separation between the NN spins ("lattice constant"). In the Umit 
T — > 0, the correlation length diverges exponentially as T — > 0: 



(2.36) 



a 2 

Equations 2.34 and 2.36 tell us that there is no long range order at any finite T and 



1 



2 4 6 8 10 
k 



Figure 2.5: Decay of the correlations in the d = 1 Ising model. 

that T = is the critical temperature of the d = 1 Ising model. Such a behaviour 
is typical for systems in the lower critical dimension: 

d>dc: Tc^O (2.37) 

d<dc: Tc = Q, (2.38) 
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therefore we conclude that d — 1 is the lower critical dimension of the Ising 
model. 

The zero-field isothermal susceptibility of the o? = 1 Ising model is, according 
to the fluctuation-dissipation theorem. 



Xt{H = 0) = ^ ^{S,S,] 



(2.39) 



In the thermodynamic limit, xt becomes: 



Xt{H = 0) 



I 
N 

I 
N 



N/2 

i k=-N/2 

N/2 

2 ^{SiSi+k) — {SiSi 

i k=0 



(2.40) 



In the thermodynamic limit (A^ — > oo), Xt{H = 0) becomes 



xAH = 0) = /? 



tanh 



J-J/^bT 



(2.41) 



The susceptibility has an essential singularity of the type exp(l/T) as T — > = 
when J > 0. When J < 0, i.e., for antiferromagnetic NN interaction, the spin 
correlation function alternates in sign and the susceptibility is finite at all T. 
The internal energy is: 

E = -dhiZldl3 = -JN tanh (3 J. (2.42) 

Notice the similarity between the internal energy and the short-range correlation 
function V{i,i + \ This is a general property of the internal energy. 
The entropy is S = {E - F) /T, 

NJ 



S{H = 0) = NkB In [2 cosh pj] - — tanh pj. 



(2.43) 



In the Umit T ^ 0, S" ^ and in the limit T ^ oo, S ^ Nks In 2, see Fig. 2.6. 
The heat capacity Ch = T{dS/dT)^ is: 

2 



CH^NkB\ 



(3J 



(2.44) 



^cosh (3 J J 

The heat capacity is finite at all temperatures (no singularity at T = 0!), see Fig. 
2.7. 
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Figure 2.6: Temperature dependence of the entropy for the zero-field d — 1 Ising 
model. 




Figure 2.7: Temperature dependence of the heat capacity for the zero-field d — 1 
Ising model. 

2.2.4 Ising Model md = 2 

The Ising model for if = in d = 2 dimensions has been solved exactly by 
Onsager. The calculation is lengthy, therefore we will first show that the Ising 
model in d — 2 has a phase transition at finite T and then quote the Onsager's 
result. 

Consider a two-dimensional square lattice with all the spins on the outer bound- 
ary "up". This simulates an arbitrarily weak external field and breaks the up-down 
symmetry of the model. To show that spontaneous ordering exists, it is enough 
to show that in the thermodynamic limit at a sufficiently small but finite tempera- 
ture there are more spins "up" than "down". Now we will estimate the fraction of 
"down" spins from the average size of a spin-down domain. 
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i- 
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-f-4- 
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-f-4-f 



-f-4 
-f-4 



Figure 2.8: A possible domain configuration of the Ising model on a square lattice. 



With the above boundary conditions, all the domain walls must be closed 
loops, see Fig. 2.8. In a wall the neighbouring spins are antiparallel, this costs 
an energy 2 J per wall segment and the probability, that a domain with its total 
wall length (perimeter) b is thermally excited, is given by the Boltzmann factor 
exp[— 2/5J6]. For large b, there are many possible domain wall configurations 
with the same wall length b. For an estimate of the number m of possible wall 
configurations, assume that on each site the wall can either be straight or make 
a kink to either side. This makes a factor of 3 for each wall segment (except for 
the last one where it is bound to close the loop). Besides, a wall can start on each 
lattice site, this makes a factor of A^. Therefore, m{b) is limited by: 



m{b) < S*"-^ N 



(2.45) 



and the average number of different thermally excited domains with perimeter b 
is: 



(2.46) 



The number of "down" spins in a domain with perimeter b is < {b/AY and 
the average concentration of all "down" spins is: 

6=4,6,8,- 



1 



E 3' 



6=4,6,8,- 
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1 / q^' \ 



48 5« ,=2,3,4 



g^(4-3g + g^) 
12(1 -g)3 ■ 



(2.47) 



where a = 2/3 J, x = 26, g = 9exp(— 2q;). For large, finite (3, q is small and 
the above ratio is < 1/2. This means that the order parameter (iV+ — N_)/N is 
greater than zero at low T and the Ising model in d = 2 (and of course also in 
d > 2) has a phase transition at finite T. 

The exact critical temperature of the (anisotropic) Ising model on a rectangular 
lattice is given by the equation: 

sinh f ^'j^ ^ sinh ( ^ = 1, 
KkBTcJ KkBTcJ 

where Jx and Jj, are the spin interaction energies in x and y directions, respec- 
tively. 



2.3 Potts and Related Models 

In the previous Section we have discussed different classical two-state models 
(which are all mapped on the Ising model). A large class of classical physical 
systems cannot be described by the Ising model. As an example consider ad- 
sorption of noble gases on surfaces, like krypton on graphite. Graphite provides 
a two-dimensional hexagonal lattice of adsorption sites. Kr is large and once it 
occupies one hexagon, it is unfavourable for another Kr atom to occupy any of its 
NN hexagons, see Fig. 2.9. At full coverage, thus, Kr forms a regular lattice with 
1/3 of hexagons covered by Kr. Since there are three entirely equivalent positions 
(states!) for Kr (denoted by A, B, and C), the system is described by a three-state 
model. The adsorbed Kr forms A,B,C domains, separated by domain walls. It is 
clear from Fig. 2.9 that the A-B and the B-A walls are not equivalent, they have 
different energies in general. 

In the Potts model there are q states available to each site. The Potts Hamilto- 
nian (with only NN interactions) is: 

n^-JJ25,„,.. (2.48) 

{id) 
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Figure 2.9: Kr on graphite. There are three equivalent adsorption sites, but Kr 
atoms are too big to occupy all of them - only one third of the sites (i.e., either 
"A", "B", or "C" sites) can be occupied. We talk about "A," "B," or "C" domains. 
The domains are separated by domain walls. The A-B domain wall has a different 
structure (and energy) than the B-A wall. 



J is the interaction energy and Sai,aj is the 5-function which is = 1 if the NN sites 
i and j are in the same state and zero otherwise, 

_ ( 1 ai^aj with (ij) e NN 
"^-^^ ~ \ otherwise ^ ^ ^ 

ai = 1, 2, • • • , g is a number which describes the state of the site i. Thus, in the 
Potts model, the interaction energy is equal to — J if the two NN sites are in the 
same state and vanishes otherwise. For q — 2, the Potts model is equivalent to the 
Ising model (5* = ±1) because each site has two possible states in both cases. The 
two states can be visualized as two points on a line (which is a one-dimensional 
object). The mapping from the Potts to the Ising model is accomplished by writ- 
ing: 

Ss.s, = (2.50) 

For q = 3 the three states of the Potts model can be visualized as corners of 
a triangle (which is a two-dimensional object) whereas for g = 4, the four states 
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belong to comers of a tetrahedron, which is a 3 -dimensional object. We conclude 
that the order parameter of the g— state Potts model is a q— component scalar. In 
d = 2, the Potts model has a continuous phase transition for g < 4 and a first-order 
transition for q> A. 

Order parameter. Let xa, xb, • • ; Xq be the fractions of cells (sites) that are 
in the states A, B, ■ ■ ■, q, respectively. Of course, xa + xb + ■ ■ ■ + Xq = 1. In 
the disordered phase, all the fractions are equal whereas in the ordered phase one 
state, say A, will have greater probability than the other, xa > {xb^xc, - ■ ■ ,Xq). 
We make a simplified approach, assuming that xb = xc = ■ ■ ■ — Xq. The order 
parameter m measures the asymmetry in occupancies, 

Xb-\ \-Xq 

m ^ XA z 

q-1 

1 -|- (g — l)m 1 — m 

Xa = , Xb = Xc ^ ■ ■ ■ ^ Xq^ 

q q 

It is easy to verify that for tha = 0: xa = xb = xc = ■ ■ ■ = Xg and for m = 1: 
Xa = ^ and xb = xc = ■ ■ ■ = Xq = 0. thus, tua is the order parameter which is 
7^ in the ordered phase and vanishes in the disordered phase. 

I would like to mention another generalization of the 2-state model, that is the 
g— state (chiral) clock model. In this model, the states are imagined as hands of 
a clock (i.e., vectors on a circle) that can take only discrete orientations, 0j = 
(27r/g)(Ti (where = 0, 1, 2, ■ ■ ■ , g — 1). The g— state chiral clock Hamiltonian is 

7i = - J ^ cos (Gi - e^- + A • fij) (2.5 1) 

A is the chirality which breaks the symmetry of the model, it produces an asym- 
metry in the interaction between the NN sites. The interaction energies between, 
say, ai — 0, aj — 1 and ai — 1, aj — are different. The two-state clock model is 
identical to the Ising model and, if A = 0, the 3-state clock model is identical to 
the 3-state Potts model. For g > 3, the clock and the Potts model are not identical. 
The two-dimensional clock and Potts models are widely used in surface physics. 
The above mentioned Kr on graphite can be described with the 3-state chiral clock 
model with A 0. 

2.4 Two-dimensional xy Model 

A very interesting situation appears when the spin is a continuous vector in a plane 
(dimensionality of the order parameter is 2) and the spatial dimensionality is also 
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d = 2. As we have seen, the Ising model in d > 2 does have a phase transition 
at finite Tc- In the xy model, however, spin excitations (spin waves) are easier 
to excite thermally than for the Ising model, they are so strong that they even 
destroy long-range order at any finite T (Mermin- Wagner theorem). On the other 
hand, there has been evidence that some kind of a transition does take place at 
a finite Tc. What is the ordering and what is the nature of this transition? This 
question was answered in the early seventies by Berezinskii and by Kosterlitz and 
Thouless. 

The Hamiltonian of the xy model is: 

n^-jJ2Si-Sj^-jJ2 cos($i - $j), (2.52) 

where the spins are classical vectors of unit length aligned in the plane. is the 
angle between the spin direction and an arbitrarily chosen axis and is a continuous 
variable. For simplicity, let us consider the square lattice and let the sum be only 
over nearest neighbours. Only slowly varying configurations ($» — <^ 1) give 
considerable contributions to the partition function, therefore we expand H: 

= ^J(fr\W^r)\\ (2.53) 

Here a is the lattice constant and we transformed the last sum over all lattice sites 
into an integral over the area. $(r) is now a scalar field. 

Possible excitations above the (perfectly ordered) ground state are spin waves. 
They are responsible for destroying long-range order at any finite T. We will not 
consider spin waves here. Another possible excitation above the ground state is 
a vortex, shown in Fig. 2.10. Any line integral along a closed path around the 
centre of the vortex will give 



W(rO-dr = 27rg => ,^^7 , (2.54) 




where q is the vorticity (q is integer; q = 1 for the vortex shown in Fig. 2.10) and 
is the polar angle to r. From the equations (2.53) and (2.54) we find the energy 
of this vortex: 



E 
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Figure 2.10: Two examples of a vortex with vorticity q = I. 

R is the size of the system. Notice that the energy of the vortex increases logarith- 
mically with the size of the system. In estimating the entropy associated with a 
vortex we realize that a vortex can be put on each lattice site therefore the number 
of configurations is about {R/aY and the entropy of one vortex is 

S ^2kB\n(^ . (2.56) 

The free energy, F = E — TS, is positive and large at low temperatures. That 
means that it is not very easy to (thermally) excite single vortices. As temperature 
increases, the entropy term causes F to decrease. F vanishes (changes sign) at a 
critical temperature: 

ksTc = ^. (2.57) 

Above Tc a large number of vortices is thermally excited. 

However, this does not mean that there are no vortices at low T. The easiest 
way to see this is to draw an analogy to a system of electrical charges in two 
dimensions. If A; is a unit vector perpendicular to the plane of the system, then 
£ = V{q(f)) X k is equal to the electric field produced by a charge of strength 
2neoq positioned at the centre of the vortex. The energy of an electric field is 

Eei^'jJ dV|£? = I / dV|V(g0) X k\', (2.58) 
whereas the energy of a system of vortices is 

Evort = i I dV|V$ X = -E,i. (2.59) 
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The energy of a system of vortices of strength qi is equivalent to the energy of a 
system of electric charges! There is almost a complete analogy between a system 
of vortices and the two-dimensional Coulomb gas. We have seen that, like the 
electrostatic energy of a single charge in two dimensions, also the energy of a 
vortex diverges logarithmically with the size of the system. However, a pair of 
opposite charges, i.e., a dipole does have a finite electrostatic energy because at 
large distances the fields of two opposite charges almost cancel. In analogy to the 
energy of two charges ±q at fi and r2 we write the total energy of a system of two 
opposite vortices: 

E2 vort = 27r In . (2.60) 

This self-energy of a vortex pair (of a dipole) is very small when the vortices are 
close together. Thus, it is difficult to excite a single vortex but it is very easy 
to excite a pair of opposite vortices. Whereas a single vortex influences the spin 
configuration on the whole lattice, the effect of a pair of opposite vortices cancels 
out on large distances from the pair so that the effect is limited to an area of the 
order \ fi — r2p. 

Since the size of the system does not appear in (2.60), but it does apper in the 
expression for the entropy, we will have a finite density of bound vortex pairs. The 
equilibrium density of vortex pairs will depend on the pair-pair (dipole-dipole in 
the Coulomb gas case) interaction which we neglected here. 

Now we have the following picture. At low but > temperature (T < Tc) 
there will be a finite density of low-energy vortex pairs ("dipoles") of zero total 
vorticity. At high temperatures (T > Tc), on the other hand, there will be a large 
concentration of single vortices. The critical temperature, we estimated above, is 
thus associated with unbinding (dissociation) of vortex pairs. 

Very interesting is the behaviour of the spin correlation function {Si ■ Sj). 
We already mentioned that spin-wave excitations break the long-range order. If 
we neglect the spin wave excitations (which have nothing to do with the phase 
transition, they do not behave in a critical way at any finite temperature), then: 

r /> for T < Tc ,^ 

J^^ • ^^-^ = 1 for T>Tc. ^^-^^^ 

This follows from the fact that bound vortices don't influence the spin orientations 
at large distances whereas they do if they are free. 

In a more realistic approach, one must take into account interaction between 
the vortex pairs, which was neglected in (2.60). Surrounding vortices (dipoles!) 
will screen the (electric) field and will thus reduce the self-energy of a vortex. 
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A sophisticated real-space renormalization group treatment of the two-dimensional 
xy model [Kosterlitz 1974] shows that the correlation length ^ diverges above Tc 
as: 

{ oc exp {h/t^'^) b ^ 1.5 (2.62) 

where t = (T — Tc)/Tc- Below Tc, however, the correlation length is always 
infinite and the correlations decay algebraically (like in a critical point)! Therefore 
we can say that we have a line of critical points extending from zero to Tc- The 
phase transition is characterized by a change from algebraic (power-law) decay of 
correlations below Tc to an exponential decay above Tc- 

Analogous to vortex unbinding transition in the xy model is the unbinding 
of opposite charges of a dipole in the two-dimensional Coulomb gas. At low 
temperatures, the Coulomb gas is composed of neutral dipoles, bound pairs, and 
is not conducting. At high temperature it is composed of free -|- and — charges 
which are conducting. 

Another analogy to the xy model are the dislocation lines in a (three-dimensional) 
crystal. At low temperatures, pairs of dislocation lines with opposite Burgers vec- 
tors can exist. At Tc, the dislocations dissociate. This is the basis of the dislo- 
cation theory of melting. In this latter example, stiffness can play the role of the 
order parameter, it is 7^ below the melting temperature and vanishes in the Uquid 
phase. 

Let me stress again that the described behaviour of Coulomb gases or xy sys- 
tems is specific only for two-dimensions\ 

2.5 Spin Glasses 

The (vitreous) glasses are amorphous materials, their lattice is not crystalline but is 
strongly disordered. Topological defects destroy the long-range order in the crys- 
talline structure whereas the short-range order is still preserved to a high degree. 
Example: topological defects on a square lattice would be triangles or pentagons. 
Of course, the defects cause also a deformation of the squares. 

Spin glasses are the magnetic (spin) analogue of the (vitreous) glasses. In 
magnetic systems, the exchange interaction is often an oscillating function of dis- 
tance between the spins. If such a system is diluted with non-magnetic ions, the 
magnetic atoms are randomly positioned on the lattice and thus experience ran- 
dom exchange interactions with other, neighbouring atoms. Example: A\X\-x^^x, 
where x is small. 
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Another system, which behaves in a similar way, is a random mixture of a 
ferromagnetic and an antiferromagnetic materials. The exchange interaction be- 
tween nearest neighbours i and j is random, it depends on the type of the magnetic 
atoms i and j. 

A Hamiltonian which reflects the essential features of a spin glass is: 

^ = - E JiRi,j)Si -Sj-Y^H- Si. (2.63) 

J{Ri,j) is a random variable and the first sum includes also more distant pairs than 
just nearest neighbours. Random exchange interactions cause that - even at low 
temperatures - the spins are not ordered, the system has no long-range order in 
the the ground state. The "usual" long-range order parameter 

^ = ^ Y.{S^) (2.64) 

i 

vanishes for a spin glass. Yet, in the experiment, a broad maximum in the specific 
heat and a rather sharp maximum in the zero-field susceptibility, connected with 
hysteresis and remanescence, are seen. So, there must be a kind of a transition 
as the temperature is varied, we shall call this a spin-glass transition and the low- 
temperature state the spin-glass phase. 

At high temperatures, the spins behave like in a normal paramagnet, they flip 
around dynamically so that the thermodynamic average {Si) vanishes. There is no 
long-range order in the system and also the local spontaneous magnetization, i.e., 
the (time) average of the spin at a site i vanishes, (5;) = for each site. At low 
temperatures, however, the spins freeze in a disordered configuration, the system 
is in a state where {Si) ^ although the average magnetization still vanishes, 
m = 0. The spin-glass transition is thus a freezing transition. 

Something happens to the dynamics of the spins, therefore we must consider 
time-dependent spins, S{t). To distinguish between the two phases, we introduce 
the (Edwards-Anderson) order parameter: 

q = lim lim ^ Y.{Sr{to)S,{to + t)). (2.65) 

c— »oo iv — »oo iV 

q clearly vanishes in the paramagnetic phase and is > in the spin-glass phase, 
see Fig. 2.12. It measures the mean square local spontaneous magnetization. 
Notice that freezing is a gradual process, which starts at T/ where the first spins 
freeze and ends at T = where all the spins are frozen. 
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Figure 2.11: Spins in the low-temperature, spin-glass phase are frozen in a disor- 
dered configuration. 



Since, in the spin-glass phase the spins are (at least partially) frozen in a dis- 
ordered configuration, we can easily imagine that some spins are frustrated and 
that the spins do not freeze in a unique way. Freezing is history dependent, the 
configuration into which the spins freeze, depends on the conditions (e.g, on the 
magnitude of the magnetic filed) close to the freezing temperature. If we, upon 
cooling, apply an external field pointing "up", there will be more spins frozen in 
the "up" orientation and will stay in this orientation even if we switch the field 
off in the low-temperature phase. This means that the systems can fall into a 
metastable configuration from which it is not able to escape if the temperature is 
low. The free energy, thus, has many, many local minima in the configuration 
space, see Fig. 2.13. In the thermodynamic limit (A^ — > oo), some hills between 
the local minima will become infinitely high and the system will be trapped in 
such a minimum forever. In the time average, thus, the system will explore only a 
part of the total configuration space. This means that the time average of a phys- 
ical quantity is no longer equal to the ensemble average (where one explores the 
whole configuration space). A system in the spin-glass phase is thus not ergodic 
in the thermodynamic limit! Since it is ergodic in the paramagnetic phase, we say 
that at the freezing transition the ergodicity breaks down. 
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Figure 2.12: Temperature dependence of the order parameter q. Upon cooling, 
the spins start to freeze at T/; below Tj, more and more spins are frozen until, 
eventually, at T = all the spins are frozen. 

2.5.1 Some Physical Properties of Spin Glasses 

The static susceptibility, Fig. 2.14, below the freezing temperature depends strongly 
on the way the experiment is performed. xiT) is largest and roughly temeprature 
independent in a "field cooling" (FC) cycle, i.e., if the field H is applied above Tf 
and then the sample is cooled at constant H below Tf. This susceptibility is re- 
versible, it is independent of the history. In contrast, the "zero-field cooled" (ZFC) 
susceptibility is obtained if the sample is cooled to T ^ at i/ = and then the 
sample heated in the field. The difference between xfc and xzfc confirms the 
conjecture that the system can be frozen in different spin configurations (which, 
of course, also have different susceptibilities). 

The conjecture that freezing takes place at Tf is confirmed also by the ac 
susceptibility measurements, shown in Fig 2.15 for Cui_^.Mn£ (x = 0.9%) at dif- 
ferent frequencies. The lower the frequency, the more pronounced (less rounded) 
the cusp is. We see also that the position of the cusp shifts to lower temperatures 
as the frequency decreases. So, the freezing temperature should be defined as the 
maximum of x{^) in the limit as — > 0. 
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Figure 2.13: Free energy landscape in the configuration space of a spin glass at 
low temperatures. There is one global and many local minima. 
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Figure 2.14: Temperature dependence of the static susceptibility of Cui_a;Mna; for 
x = 1.08 and 2.02%. (b) and (d): cooled to low temperature at if = and then 
heated at H = 5.9 Oe; (a) and (c): cooled and/or heated in H = 5.9 Oe. [From 
Nagata et al., Phys. Rev. B 19, 1633 (1979).] 
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Figure 2.15: Temperature dependence of the frequency-dependent susceptibility 
of Cui_^Mn^ (x = 0.9%. [From C.A.M. Mulder et al., Phys. Rev. B 23, 1384 
(1982), ibK 25, 515 (1982).] 
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